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Abstract 

We consider thermal properties of the susy singlet majoron model. We compute 
the critical temperature Tc and the subsequent reheating temperature Trh for 
i?-parity breaking. Succesful baryogenesis constrains the parameter space via 
the requirements that Tc and Trh are lower than the electroweak phase transi¬ 
tion temperature. A further constraint is provided by requiring that the gauge 
singlet should be kinematically allowed to decay, in order not to have a mat¬ 
ter dominated universe at the time of nucleosynthesis. We have made a detailed 
study of the parameter space and find an upper limit for the susy breaking scalar 
mass mo < 750 (900) GeV if mgin in n = 100 (1000) GeV, which is valid except 
for certain special values of the singlet sector parameters. 
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In the minimal supersymmetric standard model (MSSM), gauge invariance and 
supersymmetry allow terms in the superpotential which violate baryon(i?) or lepton(L) 
number. The absence of these dangerous terms is ensured by assuming a conserved 
discrete symmetry, called i?-parity, Rp = where S is the spin of the 

particle |^. Because of its connection to lepton number, it is interesting to consider 
situations where /2-parity is slightly broken, either explicitly p or spontaneously |^. 
Spontaneuous breaking of lepton number is possible in MSSM without introducing 
additional fields, since the scalar partner of the neutrino may acquire a non-vanishing 
vacuum expectation value 0, giving rise to a majoron ii, which in this case is mainly 
the supersymmetric partner of the neutrino. Since such a majoron should be detected 
in decay, the measurements of the width rule out this possibility, unless one 
assumes that there is small explicit L violation, in addition to the spontaneous one, 
in the MSSM [0. On the other hand it is possible to break i2-parity spontaneously 
in the MSSM with the inclusion of additional singlet superhelds P], such that the 
resulting majoron, which is now dominantly a singlet under the gauge group, is not 
in conflict with the measurements of the Z^ width. A particularly attractive scheme 
which implements such a spontaneous violation of //-parity is the supersymmetric 
version |P, of the singlet majoron model. This model is the simplest extension which 


incorporates successfully the spontaneous L-(and /2-parity-)violation in the MSSM. 
In this model there are, besides the chiral superhelds of the MSSM, right-handed 
neutrino chiral superhelds Ni (one for each generation) and an additional gauge singlet 
superheld <I>, having two units of L. The superpotential of the susy singlet majoron 
model, which is invariant under gauge symmetry and L is, in the standard notation. 


W = h^QU^H2 + h^QD^Hi + h^LE^Hi + E'LNH2 + /X//1//2 + XNN^ , 


( 1 ) 


where we have suppressed the generation indices. Eq. (|^ contains the usual terms of 
MSSM together with the Yukawa interactions for the right-handed neutrinos Ni and 
an interaction term for the gauge singlet <I>. It has been shown, through an analysis of 
the renormalization group equations (RGE), that for a wide range of parameters one 
can obtain radiative breaking of /2-parity in this model [0, 101. It has also been argued 


that the L-(and /2-parity-)-violating transition may take place after the electroweak 
phase transition. In that case the sphaleron-induced //-violating transitions are frozen 
out, so that any pre-existing baryon asymmetry is unaffected by this phase transition. 

In this paper we extend and improve the analysis of ref. |P] and study the L- 
breaking cosmological phase transition of the susy singlet majoron model in detail. In 
particular, we compute the critical temperature Tc and the reheating temperature Trh 
and compare them with the critical temperature of the electroweak phase transition. 
We also point out that the gauge singlet sector should decay fast enough in order not 
to clash with primordial nucleosynthesis. All the constraints, when taken together. 









imply that succesful baryogenesis in the susy singlet Majoron model requires an upper 
limit on the susy breaking scalar mass niQ. 

To start, we recall the tree level potential for the susy singlet majoron model [y, pi| . 
It consists of three pieces, and can be written along the neutral directions as 

Hi) + $) + N„ <h. Hi, H^) . (2) 


Here, Vh is the Higgs potential of the MSSM. We assume that the couplings \ij are 
real and work in the basis where they are diagonal, in which case we can write the 
potential {\ij = \i5ij) 


where and m| are soft susy breaking masses, and Ai are the trilinear couplings. 
Finally, if we assume that the coupling constants h'^ are small, and retain only the 
leading terms in h'^, then we may write the third piece of the potential Eq. (^) as 
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(4) 

(5) 


where m?. are the soft susy breaking masses, are the trilinear couplings and g and 
g' are the SU(2) and U(l) gauge couplings, respectively. Assuming, as in the case of 
MSSM, that all susy breaking masses and trilinear couplings are equal to a universal 
mass ruo ~ 10^ — 10^ GeV and a universal coupling A, respectively, at some GUT 
scale Mu ~ 10^® GeV, the values of the parameters in the scalar potential can be 
obtained by solving the appropriate RGEs. The effect of running of these parameters 
is to drive SU(2)xU(l) breaking through Vh{HI, H^) with (Hi) = vi = ucos/3 and 
(Hi) = V 2 = usin/3, where v = vf + = (174 GeV)^. Furthermore, for a wide 

range of parameters the nontrivial global minimum of <F) is realized in such a 

manner as to break the global lepton number and i?-parity: 


( 4 >) = 0 , {Ni)=yi, 


( 6 ) 


where typically (j) ^ i/i ^ mo- Then, nonzero vevs are induced for the sneutrinos 
through the h^-coupling which connects the ordinary doublet Higgs and lepton sectors 
to the singlet sector. If ~ mo, 0 ~ |/j ~ /i ~ mo, Ui ~ U 2 ~ and A^^'> ~ mo, 
then (z/j) ~ Mw- 

The vevs 0 and Ui are determined from the minimum of the potential We 

assume that A3 S> Ai, A2 (a minimum with A3 = Ai = A2 is not favoured by the 






solution of RGEs R OT). The minimum at scale Q is defined by 
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, 2/2 = 1/1 = 0 , 


0 = - 3 ^ 30 ;^ + 2(m4 - + Al)x - 4y43m^ , ^ 43 / 0 ; > 1 . 


(7) 

( 8 ) 


Lepton number is broken spontaneously if Eq. (|^) yields the absolute minimum, which 
is obtained when 

^3 — < x < ^3 + , m| > 0 . (9) 

Assuming universal boundary conditions at the GUT scale Mu, the parameters of the 
potential Vn<^ at the low scale Q can be written as 1,0 


m\{Q) = -mo[2 + (3 - A^)K‘^ + A^K^] , 
5 

"^^1,2 = "^0 , = 2 m|(Q) - ml , 

Ai,2(Q) = AKi , A3(Q) = AK^ , 
Ai,2(Q) = 0 , A3(Q) = \qK , 


where 


K = 


1 _L ^ \ 2 , ( 
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( 10 ) 

( 11 ) 

( 12 ) 

( 13 ) 

( 14 ) 


From the scalar potential, Eq. (Q), one can derive the mass squared matrix for the 
scalar bosons 0 in a straightforward manner. In this paper we shall only consider 
the situation when A 3 is the largest of the non-zero values of Aj, so that Eq. ( 0 ) 
holds. In the limit ^ 0 the mass squared matrix is of a block diagonal form 
= diag{M]j, M?, where Mfj is effectively the 2x2 mass squared matrix 

of the scalar Higgs bosons in MSSM M| is the 3x3 mass squared matrix of 
left-handed sneutrinos, and 


^N3<S> 


8A||/30 — 2A^\^y^ 


. 8A3I/30 — 2A3A3I/3 


AXlvl 


(15) 


To study the L-breaking phase transition, we shall focus on the (A^s, <I>)-sector only. 
At high temperatures and $ are brought into equilibrium by decays, inverse decays 
and ordinary scattering processes. The thermally averaged rate for a given process is, 
neglecting hnal state blocking, given by 


r = /n‘ini( 2 U‘‘W--P/) n 4 |m| 


initial 


(16) 

















where dllk = d^Pk/{‘^'^)^‘^Ek, fi is the momentum distribution of the particles in 
the initial state, and n* is their number density. In Eqs. (17) and (18) a constant 
matrix element is assumed. These rates are to be compared with the Hubble rate 
H ~ 23T^/Mp;. From Eq. m one easily hnds that for A not too small, N 3 and $ are 
in equilibrium already at temperatures much above 1 TeV, because of the reactions 
TVgcJ) —iVgcl) and d* —*• N 3 N 3 . The helds of the MSSM can also be assumed to be in 
equilibrium already at high temperatures. 

On the other hand, it is not obvious that the {N 3 , <h)-system is in thermal contact 
with H 2 and 9. This is because the coupling is only through h^, which is small. 
Here the most important reaction is the decay —>■ 9iNj, which induces chemical 
equilibrium at temperatures 

T< 194|-^hr//^ TeV. (19) 

If /i ~ 100 GeV, there is no thermal contact at the electroweak phase transition if 
h'^ < 10“'^. This would mean that the (V 3 , <h)-system and the MSSM would experi¬ 
ence different temperatures. In what follows we tacitly assume that this is not the 
case. For light d*, N 3 and H 2 also the scattering process is impor¬ 

tant. For instance, at m/T ~ 0.1 the particles are kept in equilibrium at temperatures 
T < lO^^Agh*^^. However, for increasing particle masses the thermalization tempera¬ 
ture decreases rapidly. 

At high temperatures the plasma masses of and <I> are given by 0 

M|(r,Q)=m|(Q) + iA=(Q)rU Ml(T,Q)^m\(Q) + \^(Q)T\ ( 20 ) 

where the RGE masses m?{Q) are given by Eq. ( 0 ). and we have neglected terms of 
order 0{h'^). Here Q is the renormalization point, which in principle is independent of 
T. However, as in thermal bath the average momenta of the particles are peaked about 
3T, it makes sense to choose the renormalization point to depend on temperature, so 
that Q ~ T. Setting Q ~ rUgusy — tuq as usual would yield masses which would differ 
by terms of the order of 0{X)] for our purposes this difference is inessential. With this 
choice, we may use Eq. (^0[) to search for the critical temperature Tc, at which the 
high T minimum (V 3 , d)) = (0, 0) becomes unstable. 

Note that the expressions Eq. ( 0 ) are valid if the temperature is much bigger 
than the masses circulating in the loops. There is a subtlety here in that at some 
point m‘X{Q) becomes negative, and the hnite T perturbation expansion, leading to 
Eq. (|20|), becomes unstable. This can be remedied by performing resummation in the 
graph involving loop, which effectively replaces by the plasma mass 

but does not change the potential to lowest order in A. As —> 0 for T —>• 
the plasma mass of may always be assumed to be less than T. Hence loops of 
always contribute to Eq. (|20|) . 



The resummed loop involving <h becomes Boltzmann-suppressed when > T. 
This implies that if tuq > 2T^, loops of $ should be neglected. Effectively, this means 
that Eq. (^Ul). The critical temperature for T-breaking is, 

however, in both cases dehned through M^(Tc,Tc) = 0. 

For a given A and A, it is then possible to run the RGEs to hnd Tc using Eq. (1^). 
In addition, one has to impose the conditions of Eq. O to make sure that a zero 
temperature global minimum exists for the chosen set of parameters. We have done 
this numerically, and the result, displaying the allowed region for hxed Tc/ttiq, is shown 
in Fig. 1. 

The parameter space can be constrained by noting that the physical eigenstates of 
the (iVa, <I))-sector should be unstable. Otherwise the energy density posc associated 
with coherent held oscillations about the vacuum would soon after the i?-breaking 
phase transition start to dominate the energy density of the universe, with disastrous 
consequences for e.g. primordial nucleosynthesis. Roughly, posc — m‘^v‘^{T/Tc)^ where 
m and v are a generic mass and vev, respectively. Thus, in the absence of dissipation, 
Posc would soon become larger than radiation energy density prad ~ T"^- 

The eigenstates can easily be found from Eq. (|l^ for a given value of A and A. 
It turns out that the heavier eigenstate, with mass m 2 , is predominantly <I>. Because 
$ does not couple directly to the MSSM states, it decays either via the small mixing 
with or via an virtual state. Here we shall assume that the latter is the case. 
Requiring that the process $ — N 3 N 3 —> where is the lightest neutralino, 

is not kinematically forbidden implies > m^o > 18 GeV, where the last 

hgure is the experimental lower bound on LSP p. The mass difference — m^y is a 
function of tuq, A and A, and the different contours are displayed in Fig. 1. If the gluino 
is light, rugiuino — 100 GeV, the experimental lower limit on sneutrino, > 37.1 GeV 
[|], rules out part of the parameter space as indicated in Fig. 1. 

Note that the actual decay rate depends on the ^ — N 3 mixing and the mixing of 
the gauge singlets with the MSSM sparticles, which is characterized by the unknown 
(but small) coupling h‘'. We have checked that $ — iVs mixing is indeed small in 
the physical large A region. It is conceivable, though, that for some values of h'^ the 
heaviest eigestate could decay directly to MSSM particles already before the onset of 
nucleosynthesis. In any case there always is some kinematic constraint that must be 
satished in order that the decay is possible. 

The critical temperature for L-breaking can further be constrained by baryoge- 
nesis considerations. L-vioIating interactions, if in equilibrium with the anomalous 
elecroweak B + L-violating interactions, will wash out any pre-existing baryon asym¬ 
metry. This can be avoided if i?-parity breaking couplings are very small, so that 
these interactions are never in equilibrium. In this case R-parity breaking would not 
be relevant for experiments. Another possibility, as suggested in |p, is to impose the 


condition that L-breaking takes place after electroweak phase transition, so that the 
sphaleron induced transitions have already dropped out of equilibrium. 

In the small approximation adopted in this paper, the critical temperature for 
the electroweak phase transition may be assumed to be the same as in the MSSM. Com¬ 
puting the critical temperature is, however, notoriously difficult task, and it is conceiv¬ 
able that it is dominated by non-perturbative effects like in the non-supersymmetric 
Because the transition is presumably only weakly hrst order, and the la- 


case 
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tent heat release is small, for our purposes it suffices to estimate the electroweak critical 
temperature by the spinoidal instability temperature Tq, determined by det^MfjiTo)) = 
0. The resummed one-loop mass matrix Mfj has been presented in |T^. The expres¬ 
sions are lenghty, and we do not reproduce them here. Using these, we have varied 
the parameters of the MSSM in the following ranges: 1 < tan/? < 50, the trilin- 
ear coupling associated with the stop-sector 0 < A* < 1 TeV and the Higgs mixing 
0 < /i < 1 TeV. We have explicitly studied two cases, the case of a light gluino with 
’^giuino = 100 GeV, and the case of a heavy gluino, rugiumo = 1 TeV. The resulting Tq 
is shown in Fig. 2. 

Succesful baryogenesis requires that < Tq. Moreover, we should also require 
that the reheating temperature Trh after i?-parity breaking does not exceed Tq so as 
to make sphaleron interactions operative again ||^ . The reheating temperature is given 
by 


T 

J RH — -5 

V 9 . 


( 21 ) 


where Unin is the global minimum and we have taken ~ 100. In Fig. 2 we show Trh 
for a choice of parameter values. In Fig. 3 we have taken all the constraints together 
to hnd the allowed parameter space. If H > 2, we hnd the upper limits 

mo < 750 GeV (rngi^ino = 100 GeV); uiq < 900 GeV (mgiuino = 1 TeV) , (22) 


where the hgures refer to the case tan/? = 1.5; if tan/? = 50, the limits are slightly 
relaxed, as can be seen in Fig. 3. However, the limits Eq. (^) are not completely 
general. If it happens that A <2 and A ~ 1, rriQ can be as high as 4 TeV. Note that 
for A -C 1 one gets H ~ 3 so that Eq. (^) always holds. We wish to emphasize that 
in all the cases an upper limit on tuq exists. 

To conclude, we have studied in detail the cosmological //-parity breaking phase 
transition in the susy singlet majoron model. In particular, we focussed on the con¬ 
straints imposed on the model by the requirements of a succesful baryogenesis. The 
necessary conditions are that T^ and Trh are lower than the electroweak phase transi¬ 
tion temperature Tq, which depends in a complicated way on the parameters of MSSM. 
The connection to the gauge singlet sector of the majoron model is provided by the 
common scalar mass parameter rriQ. An additional cosmological constraint is that the 







gauge singlet sector should be allowed to decay. This is a kinematical constraint that 
restricts the range of uiq. Because susy sparticle spectrum is essentially given by mo, 
plus D-term and radiatice corrections, baryogenesis requires a sparticle spectrum in 
the susy singlet majoron model which should be observable at LHC. 
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Figure 1; The allowed values for A and Aq. The solid lines indicate the parameter 
space for which R-parity is not broken at temperatures T > Tc with Tc/mo values 
shown in the figure. The dashed line indicates the parameter space allowed by the 
requirement of global minimum. The dotted lines are contours of m 2 — mi/mo, where 
mi ^2 are the physical eigenstates of <h and N^. From experimental lower limit for z>, 
the area above the starred-dashed line is not allowed if mgiuino ~ 100 GeV, /x = Rt = 0 
and tan /3=1.5. 








Figure 2: Spinoidal instability temperature Tq (solid and dashed lines) and reheating 
temperature Trh (dotted lines) as functions of mo- The solid lines correspond to 
^giuino = 1 TeV, fx = At = I TeV, and the dashed ones to rUgiuino = 100 GeV, 
jj, = At = As indicated, tan/? = 1.5 or 50. 








Figure 3: The allowed values of mo and A for Aq values 0.5, 1, 2, 5, and 100. The 
allowed range for mo, which is the area between the lines, is shown for the cases 
rrigiuino = 100 GeV, jj, = At = Q with a) tan f3 = 1.5, b) tan f3 = 50, and for rUgiuino = 
1 TeV, 11 = At = 1 TeV with c) tan /? = 1.5 and d) tan /? = 50. 





















